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Abstract 

The recently reportedBH coexistence of an oscillatory magnetic order with 
the wave vector Q = 0.241 A^^ and superconductivity in TmNi2B2C is an- 
alyzed theoretically. It is shown that the oscillatory magnetic order and su- 
perconductivity interact predominantly via the exchange interaction between 
localized moments (LM's) and conduction electrons, while the electromag- 
netic interaction between them is negligible. In the coexistence phase of the 
clean TmNi2B2C the quasiparticle spectrum should have a line of zeros at the 
Fermi surface, giving rise to the power law behavior of thermodynamic and 
transport properties. Two scenarios of the origin of the oscillatory magnetic 
order in TmNi2B2C are analyzed: a) due to superconductivity and b) inde- 
pendently on superconductivity. Experiments in magnetic field are proposed 
in order to choose between them, 
e-mail : kulic@ audr ey. mpi-st uttgart . mpg.de 



1 



Typeset using REVT^ 



The problem of the coexistence of magnetic order and superconductivity is a long- 
standing one and GinzburgP was the first to note the antagonistic character of these phe- 
nomenon. Further impetus in this field came after the discovery of the ternary rare earth 
(RE) compounds (RE)Rh4B4 and (RE)Mo6X8 (X=S,Se), see Ref. ^. In many of these 
compounds both ferromagnetic (F) and antiferromagnetic (AF) ordering coexist with su- 
perconductivity (S), see Refs. It turned out that the coexistence of S and AF ordering 
was realized in many of these compoundsB usually up to T = 0, while S and modified F 
ordering coexisted in ErRh4B4, HoMogSs and HoMogSes only. The reason for this was the 
antagonistic characters of these orderings. A theory has been developed and the phase di- 
agram was given in Refs. ^P,^ where the possibility of the coexistence of S and spiral or 
domain-like magnetic order has been elaborated quantitatively by including the exchange 
(EX) and electromagnetic (EM) interaction of conduction electrons and localized magnetic 
moments (LM's). It has been also demonstrated that the theory based on the EM interaction 
onlyS can not describe the coexistence problem in real systems. Note, some heavy fermions 
UPt3,URu2Si2 etc. show a coexistence of the AF and S orderings. Recently, it has been 
also found experimentally the coexistence of nuclear magnetism and superconductivity^ in 
Auln2, which was theoretically analyzed in Ref. |T2[ 

However, recent discovery of superconductivity in the quaternary intermetallic com- 
pounds (RE)Ni2B2C (RE=Sc,Y,Lu,Tm,Er,Ho and Th) has received appreciable attention, 
because of the relatively high superconducting transition temperature - 16.6 K in Lu. Band 
structure calculationslll show that the electronic spectrum is three-dimensional. Because of 
the spatial isolation of magnetic ions there is a possibility for the coexistence of a magnetic 
(M) order and S in (Ho,Er,Dy)Ni2B2C with (T^; Tm) = (8,11,6.5 K;6,7,10.5 K) respectively^. 
Tm is the ( ant iferro) magnetic transition temperature. These compounds are characterized 
by the ferromagnetic alignment in each layer, with the magnetic moments of two consecu- 
tive layers aligned in opposite directions. Band structure calculations^! of the nonmagnetic 
LuNi2B2C compound show that the conduction electron density on the RE ions is small 
(similarly as in (RE)Rh4B4 and (RE)Mo6X8 ) giving rise to a relatively small exchange 



energy. 

The subject of this paper is the theoretical analysis of the coexistence problem in 
TmNi2B2C, which is superconducting below Tc ~ 11 K with an oscillatory magnetic or- 
dering of the Tm moments below Tm ~ 1.5 K with persisting coexistence up to T = 0. 
The magnetic structure is incommensurateB with the Tm moments along the c-axis and 
with a sinusoidal modulation of their magnitudes along the (110) direction. This sinusoidal 
order is characterized by the wave vector Q = 0.241 A^^. One should stress the following 
facts: (1) TmNi2B2C is unique in the (RE)Ni2B2C family, which shows a modulation of the 
magnetic order in the (110) direction and the alignment of the moments along the c-axis; 
(2) the wave vector Q is neither large nor small {C,q^, ^ Q ^ kp), where ^q, A, kp are 
the superconducting coherence length, magnetic penetration depth and Fermi momentum 
respectively. 

In what follows the theory of magnetic superconductors - the MS theoryiH, is applied to 
TmNi2B2C and it will be shown that the competition between S and the oscillatory M order 
is predominantly due to the EX interaction, while the EM one is negligible, not only in this 
compound but also in the whole (RE)Ni2B2C family. An analysis of effects of the magnetic 
field allows us to discern between two possible scenarios for the origin of the oscillatory 
magnetic order. 

The MS theory!! considers all important interactions between LM's and conduction 
electrons: (1) via the direct EX interaction; (2) via the induced magnetic field B(r) = 
curlA(r) - the EM interaction, which is due to the dipolar magnetic field Bm(r) = 47rM(r). 
The general Hamiltonian of the (RE)Ni2B2C compounds is given by 

H = I d^r{^p\r)e{p--^A)^|Jir) + [Air)^^ ir)tcry^\r) + c.c] + + H^^^ 

+ 5: /(r - r,)i;\r)aig - l)J,^(r) + ^I^lM^} + ^[-B(r,)^?e/XBJ. + ^cf(J.)]. (1) 

i i 

Here, e(p— ^A), A(r), A, /(r), V, a, Jj and g are the quasiparticle energy, the supercon- 
ducting order parameter, the vector potential, the exchange integral, the electron-phonon 
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coupling constant, Pauli matrices, the total angular moment and the Lande factor respec- 
tively. The first three terms in Eq. (1) describe the superconducting mean-field Hamiltonian 
in the magnetic field B(r) = curlA(r), while the term Himp describes the electron scattering 
on nonmagnetic impurities. The term proportional to a{g — l)Jj describes the direct EX 
interaction between electrons and LM's, while HcFi^i) is responsible for the crystal field 
effects and magnetic anisotropy. Based on Eq. (1) and by using Eilenberger equations one 
can find the free-energy functional of the coexistence phase in terms of order parameters - 
see below. 

A. Characteristic parameters of TmNi2B2C 

The critical temperature of the transition into the oscillatory magnetic state Tm ~ 
1.5 K is small compared to the superconducting critical temperature T^. ^ 11.5 K and 
it is of the order of the exchange energy Q^x - see below. From = N{Qi)h1^ we can 
estimate the exchange interaction between electrons and LM's, which is characterized by 
/lex = -^(0)(fi' — l)nJ , where iV(0) is the electronic density of states at the Fermi level (per 
LM), n is the concentration of LM's. In absence of data on A(0) and Fermi velocity vp in 
TmNi2B2C we use the band structure valued! for LuNi2B2C, where A(0) ^ 2.4 states/eV- Lu 
atom and vpx = vpy — (2 — 3) x 10'^ cm/sec. This procedure is justified because the /-levels 
of Lu and Tm ions are weakly coupled to the conduction electrons. The decrease of Tc in 
the (RE)Ni2B2C family is scaled0 by de Gennes factor [g — 1)^ which allows to estimate 
Gex = N{0)hl^ and h^x- Namely, the Abrikosov-Gorkov formula dT^/dx ~ — '7r^Bex/2 for 
the decrease of T, in Lui_^Tm^Ni2B2C , i.e. (Arc/Ax)Lu-T„, ~ T^"" - Tj"^ ^ 5 K~ bQ^x 
gives Qex ~ 1 K and hex ~ 60 K. 

The long-range part of the EM dipole-dipole interaction between LM's is characterized 
by ©em = 27™/i^ , /i = gfisJ- The neutron diffraction measurement Jil in TmNi2B2C give 
fi ~ 5fiB, while from the crystallographic structure follows n ^ 2 ■ 10^^ cm~^, which gives 
Qem ~ 2 K. Note that Qem ~ ©ex ~ Tm. From Aq ~ 1.76 Tc one obtains ^ 250 A, while 
from magnetization measurements near Hc2 and from the slope of Hc2 near Tc (see Ref. 
it follows K = A/f ^ 7 and ^ ^ 110 A, where A ^ 0.62Al(^o/0^^^ and ^ ^ 0.85(^o0^^^- 

4 



This gives the mean-free path / ^ 50 Aand the London penetration depth ~ 500 A. 
One can say that the samples studied by Cho et a/.0 were in dirty hmit where also holds: 
{K^r/hf < 1 and (Ql)-^ < 1. . 

B. Free-energy functional of the coexistence phase: 

Since Q <^ kp the problem of interplay between S and M is treated!'! using the 
Eilenberger equations for the normal gi^{\,H) and anomalous /;^(v,R) electronic Green's 
function. They describe the motion of electrons in the EX field hex(R-) = chq sin Qz 
{hq = hexSq, Sq = |(J)|/J) - the EX interaction, and in the dipolar magnetic field 
B(r) = curl A (r) - the electromagnetic EM interaction. We present only some necessary 
results for the free-energy (per LM) F{A, Sq, Q} = Fs{A} + Fm{Sq} + Fj„t{A, Sq}, where 

F,{A} = -iiV(0)AMn^. (2) 

A is the 5* order parameter and Aq is the S order parameter in equilibrium and in absence 
of magnetism. The magnetic part Fm in the mean-field approach is given by 

Fm{Sq} = - E{[0O + QexiUQ) - m Sq,± P + I Sq,,| 

Q 

+ Ge™ I Sq,|| P +Di\ Sx,Q P + I Sy^Q H} + Fo{Sq} (3) 

Here, D > and Sq,_l, Sqj| are transverse and longitudinal (w.r.t. Q) components of Sq 
respectively, while Fq{Sq} is the isotropic part (entropy term) of the functional for isolated 
ions. ©0 = ©ez + 0em/3 + ©g^, + 6g„ characterizes the contribution of all mechanisms (long 
(Bex, ©em) " and short (0g^, 9g^)-range parts of the exchange and dipole energies) to the 
ground-state energyill. We assume that in the normal state the dipole-dipole interaction 
of LM moments leads to ferromagnetic ordering, but exchange interaction may result in 
ferromagnetic or oscillatory ordering depending on Fermi surface structure. The electronic 
susceptibility Xe{Q) = ^(O)Xe(Q) is still unknown, although it can be calculated by knowing 
the band structure of TmNi2B2C. Two scenarios for Xe{Q) will be presented below. The 
interaction part, Fint{A, Sq} = F^^^ + F^^, of the free-energy contains the EX and EM 
contributions respectively. 
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For samples TmNi2B2C studied by Cho et a/.0 the dirty limit {Q I ^o) is realized, 
as well as IhexSg/vp)^, {elAq/cY <^ 1, which allows us to find Pint 

The first term on the right hand side of Eq. (4) describes F^^^ and the second one F^^ 
respectively. Eq. (5) is derived by assuming that: (a) r^A > 1 what is indeed fulfilled 
in TmNi2B2C, where i^mA ^6 — 7; (b) the Fermi surface is isotropic - fulfilled also in 
TmNi2B2C. Note, the expression for Fj„f{A, 5*^} in the clean limit can be found in Ref. ^. 

Already on this level we can estimate the relative contribution of the EX and EM terms 
in the interaction of superconducting and magnetic subsystems in TmNi2B2C. For parame- 
ters extracted from experiments - see A, one gets r = {F[^/F[^^) ~ Qem/Qexi^iQ)^ ~ 10^'^. 
This important result means that in TmNi2B2C the EM interaction makes a negligible con- 
tribution to Fint and as a consequence the competition of S and the M order in TmNi2B2C 
is exclusively due to the EX interaction . The similar situation is realized in the whole 
family of (RE)Ni2B2C compounds where r < 10~^. Moreover, when TmNi2B2C is placed 
in external magnetic field the EX interaction plays a decisive role- see below. This means 
that the approach which is based on the EM interaction only (see Ref. |18| and references 
therein) is inadequate in explaining properties of the (RE)Ni2B2C family. However, the EM 
interaction, although much less detrimental for superconductivity than the EX one, makes 
the magnetic structure transverse, i.e. S ■ Q = due to the 9em | Sqj| p term in Eq. (3). 

C. Origin of the oscillatory magnetic order 

The magnetic free-energy in the normal state Fm{Sq} depends on the electronic suscep- 
tibility Xe{Q) = XeiQ) I glf^B^ contains magnetic anisotropy {D > 0) and the single 
ion term Fq{Sq} - see Eq. (3). These quantities determine magnetic structure in TmNi2B2C 
in absence of superconductivity. At present both are unknown and therefore in what follows 
we analyze two possible scenarios for the origin of the oscillatory magnetic structure, which 
depends on the form of Xe{Q) in TmNi2B2C: 

1. Ferromagnetic (F) scenario In this scenario Xe{Q) reaches maximum at Q = , i.e. 
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the ferromagnetic order would be realized in the normal conduction state of TmNi2B2C 
below some temperature - see Fig. la. However, in the S state it is transformed into an 
oscillating magnetic structure with the wave vector Q ^ kp - see below, and because 
a^Q^ <^ 1 (magnetic length a < kp^ ) one has XeiQ) ~ ^(0)(1 — a^Q"^)- Replacing this 
Xe{Q) in Eq-(3) and by minimizing the free-energy F with respect to Q one gets the sinusoidal 
magnetic structure at T very neat Tm (when hexSq ^ A and higher order terms in | Sq p 
are negligible ) with Qm = (7r6ea;/40oa^^o)^''^- In the presence of the magnetic anisotropy 
and by lowering temperature | Sq | grows and higher order terms in | Sq p (described by 
Fo{Sq} in Eq. (3)) become important giving rise to higher harmonics 3(5,5(5, etc.. As a 
result the striped transverse one-dimensional domain structure (Sq ■ Q = 0, Sq || 2;- axis) is 
formedliji with the magnetic energy (per LM) 

Fm = Fo{Sl} - Q^xSl + r/(^J, T)|, (5) 

where r] is the domain wall energy given hy 7] = kp^e^Sl- Here, 6^ ^ 0.6(eoI?)^/^ for 
D < Qq but {D/Qof/^ > 0.25(A;^^o)"^/^ while 6^ ^ O.30o for D > Q. This phase is 
in further called the DS-phase. In the DS-phase and at T <^ Tm the wave vector of the 
structure is given by Qds ~ '^iQex/^wkp^oY^'^- Since Q"^^^ = 0.241 A^^ and by knowing kp 
- for instance kp ^ 1 A'^, one obtains reasonable value for B^; ~ (0.1 — 0.2) K. These results 
mean that in the F-scenario the transformation from sinusoidal to the domain-like structure 
takes place around Tm, with small changes from Qm to Qds, where superconductivity and 
the domain-like magnetic structure (the DS phase) coexist. Moreover, for the given set of 
parameters in TmNi2B2C - see A, one gets that at T = one has F^s—Fm ~ — 0.3A^(0)A^/2, 
where F^s is the free-energy of the DS-phase. This means that the F-scenario for the origin 
of the oscillatory magnetic order in TmNi2B2C predicts that superconductivity coexists 
with the domain-like magnetic order up to T = 0. We pay attention that the latter result is 
independent of the scenario (F or O - see below) and it is in accordance with the experimental 
finding^ in TmNi2B2C, where S and the oscillatory magnetic order coexist up to T = 0. 
/J. 0-scenario - In this scenario it is assumed that the oscillatory magnetic order with the 
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wave vector Q would be realized in absence of superconductivity, i.e. Xe{Q) is peaked at Qq - 
see Fig. lb. At lower temperatures the magnetic anisotropy and the single ion term Fq{S'q} 
transform the structure into a domain-like one. In this case Fj„t{A, Sq, Q} is also given by 
Eq. (4), where the wave- vector Q should be considered fixed (by experiment). We point out 
that the domain-like magnetic structure in the 0-scenario is a property of the normal state 
and not of the superconducting one - see Fig.lb. Note, the ratio r ^ 1, i.e. it is small in 
both scenarios. Because in the 0-scenario one has also F^s — Fm ~ — 0.3A^(0)A^/2 then 
the domain-like magnetic structure and superconductivity coexist also up to T = 0. 



X(Q) 



X„(Q) 


/Z.(Q) 






Q 



Z(Q) 



Qm 



5C„(Q) 










Q 



(a) 



(b) 



FIG. 1. Schematic shape of the spin susceptibility of conduction electrons in the normal 
- Xniff) and superconducting - state in the case of: (a) the F-scenario, where the 

peak in Xs(Q) is due to superconductivity; (b) the 0-scenario, where the peak in Xsiff) is 
independent on superconductivity. 

Note, recent neutron scattering measurements showi clearly the presence of the third 
harmonic (SQ), at T < 1 K, with the intensity /3Q ^ 0.03/q what tells us that the magnetic 
structure in TmNi2B2C is domain-like. The smallness of J3Q can be due to the presence of 
defects in the sample which are always detrimental for a domain structureilli. The F- and 
0-scenario can not be resolved by this type of measurement. For that we need to study the 
system in magnetic field. 
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D. Gapless superconductivity 

In clean superconductors the oscillatory magnetic order can give rise to the gapless 
quasiparticle spectrunJii'i if hex > A, what is just the case in TmNi2B2C where hex > 60 K 
and A < 20. The gapless region on the Fermi surface given by the condition Q ■ = 0. In 
the DS-ph&se at temperatures where hexSq > A the density of states is given by 

M,iE)^NiO)^ElJA, E«A. (6) 

By measuring tunneling conductance, where cr(y) ~ NgiV), one could test this prediction 
which is a consequence of the EX interaction. 

E. Effects in magnetic field 

Measurements in magnetic field can discern between F- and DS-scenarios. In the F- 
scenario the critical field H^^ for the first order F-DS transition (at T <^ Tm) is obtained 
from the condition -HM{0) = -N{Q)Al/2, i.e. H^^ ^ 200 G for M(0) = 10^ G (M(0) = 
nfi is the saturation magnetization). Note, in getting H^^ we have assumed that the field is 
oriented along the c-easy axis while if it is along the hard axis it could be much higher, i.e. 
H^-b ^ . On the first sight such a small value of iff^'^ contradicts reports on the critical 
field0 in TmNi2B2C, where rather high critical field Hc2 ~ 1 T is found near Tm- Concerning 
this point one should stress the measurement Jil have been done: (1) on the polycrystals; 
(2) at fixed magnetic field H by lowering temperature. Because of possibility -ff^^, ^ H^^ 
the measurements on single crystals are desired, and because the transition at H^^ is of the 
first order, with a possibility for huge hysteresis, one should perform measurements at fixed 
T in increasing and decreasing field. Note, such a huge hysteresis is not expected in the 
0-scenario, where the critical magnetic field is determined by superconducting properties 
mainly and must be much larger than H^^ . 

In conclusion, we have found that the oscillatory (domain-like) magnetic order and su- 
perconductivity coexist in TmNi2B2C up to T = and that their competition is due to the 
exchange interaction between conduction electrons and localized moments (LM's), while the 
electromagnetic interaction is negligible, eventually helping in making the magnetic structure 
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transverse. The type and the origin of the magnetic structure in absence of superconductiv- 
ity can be resolved by measuring critical magnetic fields, where a huge hysteresis could favor 
the F-scenario, i.e. the ferromagnetic order in absence of superconductivity. It is predicted 
the gapless superconductivity in clean TmNi2B2C with line of zeros on the Fermi surface. 

One of the authors (M.L.K.) would like to thank Universite Bordeaux for kind hospitahty 
and O. Andersen, L. Hedin, Y. Leroyer, M. Mehring and V. Oudovenko for support. Work 
in Los Alamos is supported by the U.S. DOE. 
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